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We present a search for y-ray emission from the direction of the newly discovered dwarf galaxy 
Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between 
~ 2 — 10 GeV and is consistent with annihilation of dark matter for particle masses less than a few 
x 10 J GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and 
taking into account trials factors, we detect emission with p-value less than 9.8 x 10 -5 (> 3.7 a). An 
alternative, model-independent treatment of background reduces the significance, raising the p- value 
to 9.7 x 10~ 3 (2.3er). Even in this case, however, Reticulum II has the most significant y-ray signal 
of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred 
for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for 
annihilation. 

PACS numbers: 95.35.+d, 98.80.-k, 95.55.Ka, 98.56.Wm 


Dark matter’s non-gravitational interactions have pro¬ 
found implications for particle physics beyond the Stan¬ 
dard Model, motivating searches for high-energy photons 
produced via annihilation. The search for y-rays in dwarf 
galaxies [e.g.UHSDJ provides an alternative to searches in 
regions that enjoy superior statistics but suffer from com¬ 
plicated backgrounds (e.g. the Galactic center [2Tf[3Tj ). 
The observed stellar kinematics of dwarf galaxies imply 
gravitational potentials dominated by dark matter [SD- 
[33] , Many of these objects are nearby, are located at high 
galactic latitudes far from complicated emission regions, 
and possess no known astrophysical y-ray sources. Pre¬ 
vious studies of dwarf galaxies have found no significant 
y-ray emission, setting strong limits on the cross section 
for dark matter annihilation 0 UMl HU essu 

Using photometric data from the Dark Energy Survey 
(DES) [32], Koposov et al. j33] and The DES Collabora¬ 
tion et al. [33] have recently announced the discovery of 
several low-luminosity Milky Way satellites in the South¬ 
ern sky. Koposov et al. |43i report 9 new objects. One 
of these, Reticulum II (Retll), at a distance of 30 kpc, 
is the nearest dwarf galaxy after Segue 1 (Segl, 23 kpc) 
and Sagittarius (24 kpc). Retll is ~ 3 times more lu¬ 
minous than Segl, suggesting that its dark matter halo 
may be more massive than Segl’s and making Retll an 
attractive place to search for a dark matter annihilation 
signal. 

Reticulum II occupies a near-ideal location for y-ray 
analysis: 49.7° below the Galactic plane and far from 
known y-ray emitting sources (the closest source in the 
3rd Fermi Catalog [45] is 2.9° away). At energies above 


1 GeV the y-ray point spread function is significantly less 
than 1°, making source contamination unlikely. Of the 
nearby dwarfs, only Segl is further from known sources. 
The interstellar emission model provided by the Fermi 
collaboration shows that emission from diffuse processes 
is relatively uniform within 10° of Retll. 

We use Fermi-LAT data ,46] collected between Au¬ 
gust 8, 2008 and February 6, 2015. Using the pub¬ 
licly available Fermi Science Tools (http://fermi. 
gsfc.nasa.gov/ssc/) (version v9r33p0), we extract 
Pass 7 Reprocessed SOURCE class events within 10° of 
Retll using gtselect with zmax=100°, and find good 
time intervals with gtmktime with filter DATA_QUAL==1 
&& LAT_C0NFIG==1 and roicut=no. The PSF and ex¬ 
posure in the direction of Retll are found by run¬ 
ning gtselect with a radius of 0.5°, gtmktime with 
roicut=yes, gtltcube with default options, and gtpsf 
(with 17 log-spaced energies between 133.3 MeV and 
1.333 TeV, thetamax=10°, and ntheta=500). 

The search for annihilation is based on event weighting 
as discussed in [4T]. The search suffers minimal loss in 
sensitivity when including only events within 0.5° of a 
dwarf galaxy and with energies above 1 GeV (see Figs. 
3-5 of [H]). We adhere to these criteria in this analysis 
and define a region of interest (ROI) as a region of radius 
0.5° containing events between 1-300 GeV. Gamma-ray 
sources from the 3rd Fermi Catalog are assigned masks 
of at least 0.8° (the approximate PSF at 1 GeV). 

Figure [T] shows the energy spectrum derived from 
an ROI centered on Retll (red points). For each 
energy bin, the differential flux dF/dE is the num- 
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the ROI, with larger values of T providing evidence of a 
signal. In this approach, the most powerful weight func¬ 
tion for testing the background-only hypothesis is given 
by w(E, 9) = log[l + s(-E, 9)/b(E, 9 )], where s(E, 9) is the 
expected number (in a small dE , dO range) of events due 
to dark matter annihilation for the alternative hypothe¬ 
sis (signal) and b(E, 9) is the expected number from all 
other sources (background). 

The expected signal depends on the dark matter parti¬ 
cle properties (mass M, annihilation cross section (<rv)), 
the dark matter content of the dwarf galaxy (parame¬ 
terized here by the single quantity J [e.g. 31]), and the 
detector response (exposure e and PSF): 


s(E,0) 

dEd9 


(av)J dNf(E) 
8t tM 2 dE 


x e(i?)PSF(0|.E)27r sin(0). (1) 


FIG. 1: Energy spectrum of events detected within 0.5° of 
Retll (red points) with 68% Poisson error bars. Two back¬ 
ground estimates are shown: 1) the sum (solid black) of 
the Fermi Collaboration’s models for isotropic (dashed) and 
galactic diffuse (dot dash) emission at the location of Retll, 
and 2) the average intensity (gray triangles) within 3306 ROIs 
that lie within 10° of Retll and overlap neither known sources 
nor the ROI centered on Retll. The number of events detected 
from Retll in each energy bin is shown above the error bar. 
The number expected from the Fermi background model is 
shown below the solid black curve. 

ber of events divided by the width of the energy bin, 
the instrument exposure, and the ROI’s solid angle. 
Error bars indicate standard 68 % Poisson confidence 
intervals [e.g. S3 on the mean counts in each bin 
(5 bins per decade between 0.2 GeV and 300 GeV). 
The figure also shows two estimates of background. 
First, the solid black line represents a two-component 
background model that is derived by the Fermi col¬ 
laboration (http://fermi.gsfc.nasa.gov/ssc/data/ 
access/lat/BackgroundModels .html). It is the sum 
of the isotropic spectrum iso_source_v05.txt (dashed 
black line) and the diffuse interstellar emission model 
gll_iem_v05_revl. f it (dot dashed). The latter is aver¬ 
aged over the 1° region surrounding Retll (we confirmed 
that the curve does not change for any choice of radius 
within 5°). Second, gray triangles indicate an empirical 
estimate of background, showing the average intensity 
within 3306 ROIs that fall within 10° of Retll and do 
not overlap with any source masks, the central ROI, or 
the boundary of the 10° region (see Fig. |3j right panel). 
The two estimates of background show good agreement. 
Between 2 GeV and 10 GeV, the spectrum from Retll 
clearly rises above the expected background. 

To derive a detection significance we employ the follow¬ 
ing method (see [41] for details). Each event in the ROI is 
assigned a weight w(E, 9) based on its energy E and an¬ 
gular separation 9 from the ROI center. The test statistic 
T = w{Ei , Of) is the sum of the weights of all events in 


For annihilation into a final state /, dNf /dE is the num¬ 
ber of 7 -rays produced (per interval dE) per annihilation. 
We adopt the annihilation spectra of Cirelli et al. 3iJ] . 
which include electroweak corrections ED]. Note that the 
unknown J value is exactly degenerate with ( av ). 

We quantify the signal’s significance by calculating its 
p-value: the probability that background could generate 
events with a total weight greater than that observed for 
the ROI centered on Retll. We also quote “cr values”, 
CDF~ 1 (1 —p), using the standard normal CDF. 

First we compute significance by modeling the back¬ 
ground in the central ROI as an isotropic Poisson pro¬ 
cess. This procedure is justified by RetlFs location in a 
quiet region that is far from known sources and strong 
gradients (see Fig. [3| right panel). Specifically, we as¬ 
sume that 1 ) the number of background events within 
0.5° of Retll is a Poisson variable, 2) background events 
are distributed isotropically, and 3) their energies are in¬ 
dependent draws from a given spectrum. Under these 
assumptions the test statistic is a compound Poisson vari¬ 
ate whose PDF we can calculate for any weight function 
and any adopted background spectrum SB ■ There is 
no assumption that the PDF follows an asymptotic form 
such as x 2 . 

We consider four possible energy spectra for the back¬ 
ground b(E , 9). The first two are sums of the Fermi col¬ 
laboration’s isotropic and galactic-diffuse models, where 
the latter is averaged within either 1° or 2° of Retll. We 
refer to these spectra as ‘Diffuse 1’ (this is the same back¬ 
ground model shown in Fig. [l]) and ‘Diffuse 2’. The third 
is an empirically-derived spectrum (‘Empirical 1’) using 
events between 1° and 5° from Retll (excluding masked 
sources). Below 10 GeV, this spectrum is a kernel den¬ 
sity estimate, with each event replaced by a Gaussian 
with width 20% of its energy. Above 10 GeV we fit a 
power law with exponential cutoff. Finally, we bin the 
same events (30 bins between 0.2 GeV and 1 TeV) in or¬ 
der to construct a fourth possible background spectrum 
(‘Empirical 2’), where the intensity between bin centers 
is found by linear interpolation in log(intensity). Fig- 
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FIG. 2: Significance of the 7 -ray excess in the direction of Reticulum II as a function of dark matter particle mass. Left: Curves 
correspond to the result of the search in various channels (i.e. using different ways of weighting events) using background model 
Diffuse 1. The curve for e + e _ is similar to fi + fi~ , ZZ is similar to W + W~, and q represents u, d, c, s quarks and gluons. Right: 
Significance in the t + t~ channel for four different background models (see text). 


ure [ 2 ] shows the significance of the detected 7 -ray signal 
from Retll for various annihilation channels and for each 
background model. In every case, the significance peaks 
above 4cr, with little dependence on choice of background 
spectrum. 

However, it is important to consider a “trials factor” 
to account for the fact that we are searching for dark 
matter particles of any mass, i.e. conducting multiple 
hypothesis tests on the same data. As shown in Fig. 6 
of [H] , the search is not particularly sensitive to the par¬ 
ticle mass used in the weight function: ~ 3 trial masses 
suffice if the true mass is between 10 GeV and 1 TeV 
for the bb and t + t~ channels. Nonetheless, we quantify 
the trials factor by simulating large numbers of ROIs un¬ 
der the Diffuse 1 model. A p-value is found at each trial 
mass and the minimum of these p m is recorded for each 
simulated ROI. The “global” p-value p g i 0 bai is the frac¬ 
tion of simulated ROIs with p m less than that observed 
in Retll. Simulating ~ 30 million background ROIs, we 
find Pgiobai = 9.8 x 10 ~ 5 for bb and p g i 0 bai = 4.2 x 10 -5 
for t + t~. Note that the trials factor may have a more 
significant effect for a lighter final state (e.g. electrons). 

Following puna no, we also consider an entirely dif¬ 
ferent procedure for computing significance. Under this 
second procedure, we construct the PDF of T due to 
background by making a histogram of T values for ROIs 
distributed over the region surrounding the dwarf. This 
procedure is model-independent and automatically ac¬ 
counts for non-Poisson background processes (e.g. due 
to unresolved sources), an effect examined by several 
groups 

The left-hand panel of Fig. [3] shows the significance 
of Retll’s signal as calculated following the model- 
independent procedure. Compared with the Poisson- 


process model for background (see above), this proce¬ 
dure assigns less significance to RetH’s 7 -ray signal (in 
accord with ES SB HB). For example, when searching 
for a 25 GeV particle annihilating to t + t~, eight of 3306 
background ROIs have T-values larger than RetlFs (2.8er; 
other channels show similar reductions in significance). 

A trials factor for the model-independent approach is 
found by counting the number of background ROIs which 
have T values among the top n for at least one mass con¬ 
sidered (n is the rank of the central ROI at the most 
significant mass). For annihilation into t + t~, n = 9 
and there are 32 such ROIs, giving a global p-value of 
32/3306 = 0.0097 (2.3cr). The same global significance 
is found by computing what fraction of simulated Pois¬ 
son background ROIs have a minimum p -value less than 
8/3306. 

The application of this model-independent procedure 
to Retll reveals its fundamental limitation: a strong sig¬ 
nal necessarily implies that very few background ROIs 
have T larger than that of the object of interest. Thus, 
poor sampling of the large-T tail prevents a robust cal¬ 
culation of significance for the Retll signal. For example, 
had we used a 5° background region instead of 10°, zero 
background ROIs would have given a T value larger than 
Retll, indicating that the significance calculation breaks 
down when there are not enough “independent” back¬ 
ground regions. In any case, this procedure clearly iden¬ 
tifies RetlFs as the most tantalizing 7 -ray signal from 
any known dwarf galaxy (left-hand panel of Fig. |3j) . 

If the 7 -ray signal is interpreted as dark matter anni¬ 
hilation, we perform a simple exploration of the allowed 
particle parameter space. As shown in [41] . for the two 
parameters M and (av ), the likelihood ratio is related to 
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FIG. 3: Left: Significance of 7 -ray detection for annihilation into r + r for various masses, calculated using the model- 
independent procedure of El- Solid and dashed blue lines correspond to Retll and Segl (another attractive nearby target). 
Gray curves correspond to the collection of dwarfs used in m as well as the 8 other newly discovered DES dwarfs. Right: The 
Fermi isotropic+diffuse model intensity near Retll. The color corresponds to intensity normalized to the value in the direction 
of Retll (at an energy of 8 GeV — other energies are similar). A 0.5° ROI is shown at the center and the small dots show 
the centers of the ROIs used for the empirical background estimation. White x’s mark the locations of known 7 -ray sources. 
Green circles are the R.OIs which have a test statistic larger than that in the central ROI (when searching for a 25 GeV particle 
annihilating to t + t~). 


T : 


1 _^(data (M, (av)) + bg) f 

l°g- rM . j-r-T - = T ~ / s{E,0), (2 

L{ data bg J E0 


where the integral is the expected number of events in 
the ROI due to dark matter annihilation. We denote the 
right-hand side as A (M,(av)). Maximizing X(M,(av)) 
yields the maximum likelihood estimate M,(av). The 
difference 2A (M, (av)) — 2X(M, (av)) is distributed as a 
X 2 variable with 2 degrees of freedom [54] when M, (av) 
are the true values of the mass and cross section. There¬ 
fore, regions of (M, (av)) space where this difference is 
less than 2.3, 6.2, and 11.8 constitute 68.2%, 95.4%, and 
99.7% confidence regions. The y 2 behavior holds only for 
large sample sizes and it is not clear if that assumption is 
valid here. In particular, for annihilation into electrons 
or muons, where low masses are preferred, there are very 
few events above 1 GeV but below the dark matter mass. 

Figure [4] shows the derived constraints on the product 
J(av) for a number of representative channels. Although 
we cannot make a direct measurement of the cross sec¬ 
tion, the constraints on J(av), combined with indepen¬ 
dent knowledge of (av ), allow us to make a prediction for 
the dark matter content of Retll which must hold if the 
7 -ray emission is due to annihilating dark matter. In the 
t + t~ channel, for example, dividing the maximum like¬ 
lihood estimate of J(av) (Fig. [4]) by the 95% upper limit 
on (av) from m yields log 10 J > 19.6 ± 0.3, where the 
uncertainty reflects the 68 % confidence region shown in 
Fig. [4j For comparison, Segl has log 10 J = 19.3±0.3 35] , 



FIG. 4: An exploration of a dark matter interpretation 

of the observed 7 -ray excess for four representative anni¬ 
hilation channels. J = Jig 10 19 GeV 2 cm -5 and (av) = 
(av )-26 10 - 26 cm 3 sec -1 . The data constrain only the product 
of J (a v) since the dark matter content of Reticulum II is cur¬ 
rently unknown. Contours represent 68 %, 95%, and 99.7% 
confidence regions. Note that this figure does not quantify 
which annihilation channel is preferred by the data, i.e. which 
channel provides the best fit to the 7 -ray spectrum. 


While RetH’s 7 -ray signal is tantalizing, it would 
be premature to conclude it has a dark matter origin. 
Among alternative explanations, perhaps the most mun¬ 
dane is the possibility that an extragalactic source lies 
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in the same direction. Computing T as a continuous 
function of sky position reveals the peak T-value to oc¬ 
cur 0.083° from the optical center of Retll [43] , an off¬ 
set similar to typical localization errors for weak, high- 
energy sources in the 3rd Fermi Catalog jJS]- Thus the 
emission is consistent with originating from Retll’s loca¬ 
tion. Searching the BZCAT [55] and CRATES |56j cat¬ 
alogs reveals a CRATES quasar (J033553-543026) that 
is 0.46° from Retll. Further work must be done to de¬ 
termine whether this particular source contributes to 
the emission, though we note that flat spectrum radio 
quasars rarely have a spectral index less than 2 [57U5H]. 
Other diagnostics, such as color-color diagrams, multi¬ 
wavelength surveys, and variability searches, may even¬ 
tually reveal the presence of active galaxies behind Retll. 
These must then be considered as possible 7 -ray emitters. 
However, we emphasize that even without knowledge of 
specific background objects, the p -value derived from the 
background sampling procedure (Fig. |3j) automatically 
accounts for the probability that a chance alignment is 
causing RetH’s 7 -ray signal. 

There is also the possibility that 7 -ray emission arises 
from within Retll, albeit through conventional processes. 
One of the much-discussed astrophysical explanations for 
the apparent Galactic Center excess is millisecond pul¬ 
sars [21112S I55U55] . In the case of Retll, it is the high- 
energy behavior of the emission which disfavors a pulsar 
model, as millisecond pulsars exhibit an exponential cut¬ 
off at around 2.5 to 4 GeV [2511551I55H55] . Alternatively, 
high-energy cosmic ray production could potentially arise 
in the vicinity of young massive stars. Upcoming photo¬ 
metric and spectroscopic analysis of Retll will check this 
possibility. 

Thorough explorations of the diffuse background, the 
7 -ray events toward Retll, properties of RetH’s dark mat¬ 
ter halo, and any coincident sources will prove crucial to 
confirming or ruling out the dark matter interpretation. 
Fermi’s upcoming Pass 8 data release E3[5H] will improve 
every aspect of the instrument response, allowing for a 
more sensitive analysis of Retll and other known and 
as-yet-unknown Milky Way companions. Understanding 
the 7 -ray emission, along with the analysis of Retll as a 
galaxy embedded in a dark matter halo, may provide a 
long-sought avenue for the characterization of dark mat¬ 
ter particles. 

We note that the Fermi collaboration has simultane¬ 
ously performed an independent search for 7 -ray emission 
and reports no significant excess from any dwarf galaxy, 
including Retll m can- Nevertheless, the strongest sig¬ 
nal they find (p = 0.06), for any annihilation channel and 
mass, corresponds to a 25 GeV particle annihilating into 
t + t~ in Retll (cf. our Fig. [ 3 ]). The reason for any dis¬ 
crepancy with our result is unclear, as the Fermi analysis 
is based on unreleased data. 
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